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^^ . Abstract 

^ ■ We show how to obtain the dual of any lattice model, with inhoniogeneous local 

interactions, based on an arbitrary Abelian group in any dimension and on lattices with 
arbitrary topology. It is shown that in general the dual theory contains disorder loops on 
the generators of the cohomology group of a particular dimension. An explicit construc- 
ts^ I tion for altering the statistical sum to obtain a self-dual theory, when these obstructions 
exist, is also given. We discuss some applications of these results, particularly the ex- 
istence of non-trivial self-dual 2-dimensional ^jv theories on the torus. In addition, 

A ■ we explicitly construct the n-point functions of plaquette variables for the U{1) gauge 

^ ■ theory on the 2-dimensional g-tori. 



m 
(N 






X 



PACS codes: ll.15.Ha, 12.40Ee, 75.10.H, 02.40Re. 

Keywords: duality, topology, lattice theories, spin models, homology. 



^This work is supported in part by NSERC of Canada, NATO CRG 970561 and the Fonds zur Forderung 
der Wissenschaftlichen Forschung in Ostcrreich under project number J01185-PHY. 



It has been known for half a century [jl[ that duahty is a powerful tool for analyzing 
statistical models 0]. In recent years, these ideas have been revived and have had impact in 
the study of supersymmetric gauge theories and superstring theory. The basic idea is that, 
given a statistical, field theoretical or string model, there may exist another representation 
of the model in which strong and weak coupling limits are interchanged. Most discussions 
of duality in statistical systems apply to spaces which have trivial topology. However, there 
are many instances of physical interest, such as regularized finite temperature gauge theory, 
where a statistical system is defined on a space with non-trivial topology and where some 
version of duality could be a useful tool. In this Letter we shall show that duality can indeed 
be implemented on a topologically non-trivial space. Our results apply to models with degrees 
of freedom in Abelian groups which live on /c-cells (sites, links, plaquettes, etc.) of a lattice 
that is the triangulation of a d-dimensional smooth manifold. We restrict our attention to 
models with nearest-neighbor interactions. 

The essence of our result is that Abelian duality on a lattice with non-trivial homology 
requires the appearance of disorder defects on cohomology cocycles of the lattice in either 
the original or the dual theory. These defects are similar to the familiar t'Hooft loops which 
are used to characterize the phase structure of gauge theories. The defects can be classified 
systematically and this allows one to identify self-dual models on spaces with non-trivial 
topology. For example, given a model that is self-dual on an infinite space, one can identify 
the modification of the statistical sum which makes the model self-dual when some of the 
dimensions are compactified. We will present as an example, a self dual Ising model and a 
self-dual Zjq spin model on the 2-torus. Also, we explicitly compute the n-point functions of 
plaquette variables for the U{1) gauge theory on the 2-dimensional (/-tori. 

To present our results in their most general form, we must introduce some terminology. 
We use the language of simplicial homology ||^, Q . Consider a lattice Vt and associate to every 
fc-dimensional cell of the lattice an oriented generator c^ where i indexes the various cells of 
dimension k. These objects are used as generators of the /c-chain group, denoted by Cfc(r2, G), 






Y^ giC^^^=geCkin,G) , g, E G 



Here G is an arbitrary Abelian group with group multiplication implemented through addition 
and Nk is the number of fc-cells in the lattice Q. An element g G Ck{^, G) is called a G- valued 
fc-chain or simply a fc-chain. Clearly Cfc(f2, G) = (Bi=iG. 

Two homomorphisms, the boundary d and the coboundary S, define the chain complexes 

{G,{n,G),d) and {G,{n,G),6) where G,{n,G) = ©LoCfc(fi,G): 

^ Gd{n, G) ^ . . . ^ Gfe+i(fi, G) ^ Ck{n, G) ^ Gk-i{n, G) ^ . . . Go(fi, G) ^ 
^ G,{Q, G) ^ . . . ^ Gfe+i(fi, G) ^ G,{Q, G) ^ Gk-i{Q, G) ^ . . . Go(fi, G) ^ 

These homomorphisms are defined by their actions on the generators c^ (we display the 



dimension subscripts on dk and 5^ only when essential), 



Nk-i Nk 



+1 



(0 _ Y^ uW . M) 1 ^(i) aJO _ Y^ [^(i) . Si)] JJ) 

k+l ■ ^fc J ^fe+l 



5cr = E [4^^ : 4^ii] 4^ii , 5cW = ^ [c 



i=i i=i 

where the incidence number is given by, 

|. (j) (j) , J ±1 if the j'th {k — l)-cell is contained in the i'th /c-cell 
^'"'^ ■ ""'-1^ "jo otherwise 

The plus or minus sign reflects the relative orientation of the cells. The boundary (co- 
boundary) chains and the exact (co-exact) chains are defined as. 



Bk{n,G) -- 


= Im dk+i 


, B'{n,G) - 


= Im 4-1 


Zk{Q,G) -- 


= ker dk 


, z'^in,G) -- 


= ker 6k 



These sets inherit their group structure from the chain complex. The boundary and co- 
boundary operators are nilpotent: dd = and 66 = 0. The quotient groups, 

Hk{n, G) = Zk{Vt, G)/Bk{Q, G) , H^Q, G) = Z^Q, G)/B''{Q, G) 

are the homology and cohomology groups respectively. The elements of the homology group 
are those chains that have zero boundary and are themselves not the boundary of a higher 
dimensional chain, whilst the elements of the cohomology group are those chains that have 
zero coboundary and are themselves not the coboundary of a lower dimensional chain. The 
inner product on Ck{fi, G) is defined by its action on the generators, 

and is linear in both arguments. The coboundary operator and boundary operators are dual 
to each other with respect to this inner product, 

(5c«,cif2i) = (c«,4^,) (2) 

Using this language, the models which we consider can be stated as follows: the degrees of 
freedom take values in G and are defined on {k — l)-cells, and the Boltzmann weights are 
defined on the coboundary of a (A; — l)-chain., 

Z= E X{B,{{59^cf)) (3) 

Here, and throughout our discussion, we allow the Boltzmann weights, -Bj, to differ on each 
fc-cell. This generalization allows one to consider a large class of models, for example some 
of our results are easily generalized to random bond models. Also, it enables one to use the 
partition function as a generating function for certain correlators. Of particular interest is 
the correlator of a disorder operator and an anti-disorder operator. It is obtained from the 
partition function by shifting 6g by a chain, /, which has support on an assembly of adjacent 



{d — A;)-dimensional objects denoted by F. F is constructed by choosing two arbitrary {k — 1)- 
cells on the dual lattice and connecting them by an arbitrary path of fc-cells. It is then 
re-interpreted on the original lattice by taking the lattice-dual of the path. The resulting 
object is F and the chain / is explicitly given by the expression a J2ier ^d-k ^^ some a G G. 
Clearly, there are as many disorder operators as there are elements in the group 0. 

A familiar example of a spin system is the Ising model which, in our notation, has G = ^2 
and k = 1 so that for g G Co{Q, G), 

No 

where ii{l) and ^2(0 ^^^ the end points of the link /. The Boltzmann weight for aligned spins 

( gii(i) - 912(1) = inod 2) is e^ and for anti-alinged spins ( gi^(^i) - gi,,{i) = 1 mod 2) it is e~^ 
so that, 

Z= Yl II^^pI/^ (l-2(^i,(i) -^i2(j) mod2))} 

{g^=0,l} I 

More general models of this kind include Abelian gauge theories, spin models in arbitrary 
dimensions, plaquette models, etc. 

We shall expand the Boltzmann weight in (|^) in terms of the characters, Xnio) of the 
irreducible representations R & G*, 

B,{g) = Y: h{R)xR{g) , h{R) = ^ E xAh) B,{h) (4) 

where |G| is the order of the group. In the case of a continuous group the normalized sum over 
group elements is replaced by the Haar integration measure. Since G is Abelian, G* inherits 
an Abelian group structure, where the product (taken to be addition) is implemented via the 
tensor product of representations of G. In particular: aR + bS = (®°Li-R) (8) {^'l^iS) for 
a,h E ^ and R, S & G*. This implies the following factorization properties for the characters. 



(5) 



XR{hi + h2) = XR{hi) XR{h2) R e G* and hi, h2 e G 

Xn{a hi) = XaR{hi) = X(^i^^R{hi) a e Z 

We insert (^) into the partition function (|) to obtain, 

z = E HE b^{n) x4('5^,cl'^) 

gGCfc_i(a,G) i=l neG* 

Nk ( N 

= E E X{\h{{r.cf)) x^,.^i.^{{5g,cf))] (6) 

<?GGfe_i(n,G) reCk(^,G*) i=l^ *" ^ 

Here, the product over fc-cells has been interchanged with the sum over group representations. 
Every fc-cell has been associated with a representation r^ G G* and this information is encoded 
in the G* valued fc-chain, r = '}2i=i'^iCk ■ Applying the factorization properties (1) to the 



product of characters in (^ one finds, 

i=i i=i *= \ j=i J i=i j=i 

The representations in the subscript of x can be rewritten in a more transparent manner by 
noticing the following, 

which follows directly from linearity of the inner product and the definition of the coboundary 
operator (0). It is now possible to sum expression (|^) over g, 

E n x^sr,^^^ (9,) = n E x^,,^^^ (9,) = i^r-^ n ^ ((^^ ^i^^,), o) (s) 

geCfc_i(C,G) i=l i=l 3,eG i=l 

The last equahty follows due to the orthogonality of the characters, |G|^^ ^gec Xr{.9) Xsi9) — 
6{R, S) for R, S E G*. The partition function now depends on /c-chains with constraints, 

Nk A^fc-i 

However, the constraints simply force r to be an exact chain [dr = 0). This implies that r 
is the sum of a boundary chain and representatives of elements of the homology group under 
inclusion into the chain group, 

r = dr' + '£ha7a , r'eCk+i{n,G*) (9) 

a=l 

where {7a : a = l,...,Ak} are the generators of Hk{^,G*) = (Ba=iHk,a{^,G*) and ha G 
Hk,a{^,G*) "-^ G*. We will use h E Hk {Q,G*) to denote this inclusion, which is necessary 
so that addition in the argument of the Boltzmann weight is group multiplication in G* and 
not in Hka {^, G*). Removing the constraints and inserting this form for r into the partition 
function one finds, 

Nk 

Z = |G|A^.-.-<^..i J2 E n ^^ ((^ + dr, c«)) (10) 

heHk(n,G*) r€Ck+i(n,G*) i=l 

Here we have dropped the prime on the dummy index r', and dk = dim ker d^ which is 
required to prevent overcounting. 

The final step in the duality transformation is to re-interpret c^ as elements on the dual 
lattice: so that we make the association c^_^ ^^ c^ where c^ live on the dual lattice ^2*. 



Under this association {dc^ , q ) = {5c}_f^, c2_i) follows naturally. With this reinterpretation 
the dual partition function is given by the compact expression, 

Z=|Gr™-^-- E E n h{{h + 5r,c*}%)) (11) 

where D^_f^ = d^ and is the dimension of the kernel of 6d-k on the dual lattice. Notice that h 
is now an element in the {d — A;)'th cohomology group of the dual lattice, which is isomorphic 
to the fc'th homology group of the lattice. 

We now see that, in general, the dual theory has additional disorder terms corresponding 
to the generators of the cohomology group of the dual lattice. By comparing ( pUj) and @, it is 
clear that, in order for (|) to be self-dual, Hk{fi, G*) has to be trivial. Of course, in addition, 
the lattice must be self-dual, k must equal d — k so that the original and dual degrees of 
freedom live on the same type of cells and G must be isomorphic to G*.f\ We interpret this 
duality in the general sense where the Boltzmann factors Bi{g) for each group element g can 
be regarded as independent coupling constants and the dual Boltzmann factors bi{r) for each 
r E G* = G are the images of these constants under the duality transformation. In the normal 
sense of duality one requires a more specific manner in which B{g) transforms into 6(r), in 
particular, one requires these to have the same functional form. 




Figure 1: This diagram depicts the cycles (7j) and cocycles (a*) of dimension 1 on the torus. 

Of course, the Ising model on an infinite plane satisfies all the requirements of self-duality. 
However, suppose one compactifies both directions so that the model now lives on a torus. 
Topology is the only obstruction to that model being self-dual. To see this explicitly we write 
down the partition function for the dual model. Firstly note that the generators of H^iVt^ G) 
are the cocycles h^ = Y.iea'^ Ci and h"^ = Z^ieo-^ Ci , where the sums are over all links in a^ 
shown in figure 1. These cocycles are the duals (in the sense of the inner product) to the 

^G = G* for any finite Abelian group. An example of a group where this does not hold is U{1) where 
t/(l)* ^ Z. 



cycles 7, which generate the homology group. Using (|Tl]), the dual of the Ising model on the 
torus has partition function, 



hiM=0,l {ri=0,l} 



J2 (1-2 (r,^(,) - ri,(i) mod 2 



Z^o-iUo-2 



+J2{^ ~ 2 (ri,(i) - Ti^^i) + hi mod 2)) + ^(^1 - 2 {vi^^i) - Vi^^i) + /12 mod 2 
ZGo-l zeo-2 

Here (3* = —\ lntanh/3 is the well known dual coupling constant, also notice that a^ fl cr^ = 
thus all links are accounted for in the Boltzmann weight. The dual partition function is then 
the sum of four copies of the Ising models in which the couplings along a* are taken to be 
(/?*,/?*), (-/?*, /?*), (/?*,-/?*) and (-/?*,-/?*). The extra terms have the interpretation of 
disorder defects [^ ||. 

We have shown that the dual theory on a lattice with non-trivial cohomology contains 
extra topological terms. In order to formulate self-dual models in these cases it is necessary 
to begin with a theory that contains a subset of the cohomology generators. We shall see 
that under duality these extra terms cancel some of the topological terms that would appear 
in the dual theory. Such models have partition functions given by, 

iVfc 

Z= Y. E Y{B.{{{Sg + h\cf)) (12) 

where H^^JQ, G) is a subgroup of H''{Q, G) generated by a subset of the generators of the 
cohomology group, {a" : a G A C {1, . . . ,^4^}}. An element h G HhJfl,G) is written as, 
h = J2aeA h"'0'"' with h"" G H^'°-{Q, G) ■— > G. As before, the inclusion is necessary so that 
addition in the argument of the Boltzmann weight is group multiplication in G rather than 
mH'''''{n,G). 

Once again the first step in the duality transformation is to perform a character expansion 
of the Boltzmann weights, 

Z= Y E E \{H{r,cf))x^^^^,^[{{59 + h),cf)) 

reCk{n,G*) heH^^^in,G) 9eCfc„i(n,G)*=i " 

In the above we have encoded the representations that live on fc-cells in the /c-chain r, and 
introduced the character coefficients bi ({r,cj^ )) as in equation @. The properties of the 
characters (^ factorize the partition function. 

The sum over g was performed previously and produced a delta function forcing r to be an 
exact chain (see (^ and (^). The sum over the cohomology elements will force additional 

6 



constraints on the representations. Using the factorization properties of the characters and 
the exphcit representation h = J2a(^A h°'0'"' one finds, 

a£A aeA 

Performing the sum over H^'"^ and product over fc-cells yields, 

E n x,.^), {{h, c«)) = n E n x,.c^))(.^c^)) i^n 

= n E xv-. (nc^')(.^c^)) (^'^) = n E x(n.") (/^'^) (13) 

Now interpreting the inner product {r,a°') as a representation of H'''°'{Q,G), which can be 
done since H^'°'{Q, G) is a quotient subgroup of G, and performing the sum over h"" forces 
these representations to be the trivial ones. Then ([T3| ) reduces to a product of delta functions. 

Putting this information together we obtain the very compact expression for the partition 
function, 

a'£A r£Ck{n,G*) \«=1 / a€A 

It is possible to solve the constraints and remove the delta functions. Since r is forced to be 
exact, take it to be of the form (^). The other constraints then allows one to solve for some 
of the coefficients hb. On lattices with no boundaries or torsion there exists an isomorphism 
between cocycles and cycles induced by the inner product. The generators can then be paired 
in the following way: (7^, a") = 1 for a = 1, . . . ,Ak (the inner product is unity since one can 
always choose the canonical generators which intersect on only one cell; compare figure 1.) 
Restricting ourselves to lattices that have these properties and solving the other constraints 
leads to, 

{dr\ a") +Y.h {lb, a") = => ha = - {dr', a") = - (r', (5a") = for a G A 
6=1 

the last equality follows since {cr""} are cocycles. Thus we see that the extra constraints force 
the coefficients of the cycles " dual" to the cocycles to vanish. Inserting these constraints into 
the partition function one obtains. 



z = \Gf'"^-'^+^Y[\H''^'{n,G)\ J2 E n^^(((^^ + ^)' 






a'eA 



h<^H 



*^^(n,G*) ''GC'fc+i(f^.G'*)*=l 



In the above A denotes the compliment of the set A, so that the homology subgroup is gener- 
ated by those cycles that have no counterpart in the modified model (p!2D . Now interpreting 

7 



this partition function on the dual lattice by making the associations mentioned earlier we 
find. 



N* 



d^.-.+i-^.-.-i II \H,_kA^,G)\ J2 



E I[b^({iSr + h) 



a'£A 



hGH'L''{Q*,G*) reCa-k-i{^* ,G*) i=l 






The generators of H -j, (Q*,G*) are the cocycles on the dual lattice that are associated with 
the set of cycles {7^ : b G A} on the original lattice. In the case of the Ising model with one 
cocycle, {cr^}, added in at the start the constraints eliminate {71}, so that {72} survives. The 
dual theory then has an extra term supported on the lattice-dual to {71}, denoted by {7*^} 
however, that cocycle is {cr^}. 

In order to illustrate the physical content of the ideas in this paper, we discuss the ^n 
model on the torus with topological terms in some detail. We will denote elements of Z^ by 
n, the representations of Z^ are isomorphic to ^tv and as such we label them by the integers 
r = 0, 1, . . . , A^ — 1. The characters are given by Xr(^) = exp{irn27r / N) . 

We discuss the most general case on the torus by choosing one of the subsets of the 
generators of the first Homology group: Si = {0}, S2 = {o"^}, S3 = {a^} or S4 = {a^,a'^} 
(compare figure 1.) For a given choice of Sa define the model by the following partition 
function, 

N-l N-1 

Z(Sa, f3) = J2 n ^ (^n(0 ~ ^i2{i) modiV) J| E IT ^ (^n(0 ~ "-^2(0 + ^<t modiV 

{ni=0} i^Sa '^GSa {h,y=0} leer 

for a = 1, 2, 3 or 4. In Villain form the Boltzmann factor is given by (equal on all links) 



.(.,^Ee.p{-f(^-.r} 



(14) 



The coefficients b{r) of the character expansion are given by, 
1 ^-1 1 



h{r) 



/^ V .-^(#— )^ 



/5 



E 

m£ZS 



e 2 ^N 



(15) 



Here j3* = (27r)^A^^/5^^ is the dual coupling constant. Thus we see that b{r) has the same func- 
tional form as B{n). By applying the above discussed rules on the cancellation of the cycles 
we find the following transformation properties of Z{Sg„j3) under the duality transformation. 

No 4 



z{s,,P) = (— j E^a,.^(sb,r) 



where the transformation matrix Ma b is given by 







/O 








1 \ 


M = 


1 
' 2 






2 




2 










u 








0/ 



The matrix M obeys M^ = 1 as should be since applying the duality transformation twice 
gives the original model. A similar structure for the case of the homogeneous Ising model in a 
different approach was discussed in []^ . The appearance of cocycles in duality transformations 
was also appreciated by Rakowski and Sen |^, however they did not discuss self-duality. 

Self-dual and anti-self-dual models can be constructed using the eigenvectors of M. Since 
M^ = 1, the eigenvalues are either +1 or —1. In particular three eigenvectors have eigenvalue 
1 giving rise to the self dual models Z(s2,/3), ^(83,/?) and Z{s-i,[3) + 2Z(s4,/3). The fourth 
eigenvector with eigenvalue —1 corresponds to the anti-self-dual model — Z(si, /5) +2Z(s4, (3). 

This strategy of choosing a complete set of cycles, computing the corresponding dual 
theories and finding eigenvectors of the transformation matrix M is a method for finding 
self-dual theories which can easily be implemented for other models and in higher dimensions. 
In general, M contains a single non-zero entry in each row. Since M is non-singular (in fact 
M^ = 1), this implies that every column also has only one non-zero entry. 

As a second example, we show how one can use the inhomogeneous form of the partition 
function to compute correlators. In particular we compute some correlators of lattice gauge 
theories on 2-dimensional (yf-tori. In this case we choose A; = 2 so that the degrees of freedom 
live on links of the lattice, and leave the group G arbitrary for the moment. The partition 
function is. 

Here / and p labels the links and plaquettes of the lattice respectfully, and we have reverted 
to the multiplicative notation for group multiplication. Since our formalism allows for inho- 
mogeneous Boltzmann weights it is possible to compute the expectation value of arbitrary 
n-point functions of plaquette variables in arbitrary representations. This includes eg. the 
filled Wilson loop. We alter the weights of the Boltzmann factors for the plaquettes in the 
support r of our n-point function as follows: 

Bp{g) -^ XSp{9)Bp{g) = E bpiR)XR{9)XSp{9) = E Wi^ " Sp)xr{9) , for all p G T 

ReG* ReG" 

This introduces a plaquette variable in representation Sp for all plaquettes p in T. The dual 
partition function is rather trivial since there are no cells of dimension negative one, and 
Hq[Q, G*) = G* . Using this fact, the n-point function is given by, 

/-rr . A T.h^G*Y\riT'Hh)Y{,^rHh- s*) 

here i labels the sites of the dual lattice, V* is the lattice-dual to V and S* = Sp. In the 
particular case of G = f/(l), and choosing the heat kernel action, 6i(r) = (^/27FJ3)~^e~^ /^^, 
one can express the partition function in terms of ^-functions , 



\per 




Ehe^ exp{-i2/3)-^N^h^} 



exp 




*\2 




n f 2(3 . 1 Y^ C* 



(16) 



For the (filled) Wilson loop, F is a collection of adjacent plaquettes with all Sp = 1. The overall 
factor in ([l|) is given by exp (-(2/3A^*)-M(F) {N* - A(F))), where A{r) is the number of 
plaquettes in F. Note that this expression is symmetric under interchanging the interior and 
exterior of F. In the infinite volume limit, Nq -^ cxo, this reduces to the familiar area law (the 
^-function factor in ( [T6| ) approaches 1). Using these techniques it is possible to analyze the 
sector of the gauge theory in which there are homologically non-trivial Wilson loops. 

These few examples by no means exhaust the applications of our formalism. A systematic 
study of its implications for confining gauge theories in higher dimensions is in progress. 
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